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Spatial fluctuations of the local electric fleld induced by a constant applied electric field in composite 
media are studied analytically and numerically. It is found that the density of states for the fields 
exhibit sharp peaks and abrupt changes in the slope at certain critical points which are analogous 
to van Hove singularities in the density of states for phonons and electrons in solids. As in solids, 
these singularities are generally related to saddle and inflection points in the fleld spectra and 
are of considerable value in the characterization of the field fluctuations. The critical points are 
very prominent in dispersions with a regular, "crystal-like", structure. However, they broaden and 
eventually disappear as the disorder increases. 
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In the study of heterogeneous materials, the prepon- 
derance of work has been devoted to finding the effective 
transport, electromechanical and mechanical properties 
of the material , which amounts to knowing only the 
first moment of the local field. When composites are 
subjected to constant applied fields, the associated local 
fields exhibit strong spatial fluctuations. The analysis 
and evaluation of the distribution of the local field has 
received far less attention. Nonetheless, the distribution 
of the local field is of great fundamental and practical im- 
portance in understanding many crucial material prop- 
erties such as breakdown phenomenon |^ and the non- 
linear behavior of composites H]. Much of the work on 
field distributions has been carried out for lattice models 
using numerical and perturbation methods Re- 
cently, continuum models have been also addressed using 
numerical techniques 

In this Letter, we study the local electric field fluctu- 
ations by analyzing the density of states for the fields. 
To illustrate the procedure, we evaluate the density of 
states for three different continuum models of dielectric 
composites: the Hashin-Shtrikman (HS) construction 
periodic and random arrays of cylinders. It is found that 
the density of states for the fields exhibits sharp peaks 
and abrupt changes in the slope at certain critical points 
which are analogous to van Hove singularities in the den- 
sity of states for phonons and electrons in solids. This 
analogy is new and powerful, and places the study of 
field fiuctuations in composites on the firm foundation of 
solid-state theory. In the case of the Hashin-Shtirkman 
construction, we obtain an exact analytical expression for 
the density of states. We first describe the basic equa- 
tions and then determine the density of states for the 
aforementioned examples. 

Consider a composite material composed of (n — 1) 
isotropic inclusions with dielectric constants and vol- 
ume fractions (/i^ (i = 2, . . . , n) in a uniform reference 
matrix of dielectric constant ei with volume fraction 
Clearly, the local dielectric constant at position r 
is e(r) = X]"=i ^i-^^'Hi"); where /'^*^(r) is the character- 



istic function of phase i which has non-vanishing value 
(r) = 1 only if r lies inside the volume Vi occupied by 
phase i. Let Eo(r) denote an applied electric field. The 
local electric field E(r), and the dielectric displacement 
D(r) are related via the relation D(r) = e(r)E(r). The 
potential field u{r) is related to E by E(r) = — VM(r). 
The local fields are obtained from the solution of the 
governing relation V • D(r) = subject to appropriate 
boundary conditions. 

Since the local dielectric constant of the composite ma- 
terial is a piecewise continuous function, we can solve the 
following equivalent equations: 

V\{r) =0, re V,, (1) 

u(r)|r- =w(r)|r+, re^V^, (2) 

ei5„u(r)|r- = ej5„u(r)|r+, r e dV^, (3) 

where i — 2, . . . ,n. The index j denotes neighboring 
phases in contact with a given inclusion i. We donote 
the outward normal derivative to the interface dVi by 
dn and the interface points approached from inside or 
outside inclusion by r^, and r^, respectively. The nu- 
merous interfaces between the inclusions and matrix are, 
in general, irregular and randomly distributed in space. 
In most cases solutions of (|l|)-(^) can be found only nu- 
merically. 

In order to show the salient features of the field distri- 
bution, we first consider an analytically tractable model 
of composite media: the HS composite-cylinder construc- 
tion [^. Although the effective dielectric constant eg of 
this model is known exactly, its local field distribution 
has heretofore not been investigated. The HS two-phase 
model is made up of composite cylinders consisting of a 
core of dielectric constant £2 and radius a, surrounded 
by a concentric shell of dielectric constant ei and radius 
b. The ratio {h/af' equals the phase 2 volume fraction 
<f>2 and the composite cylinders fill all space, implying 
that there is a distribution in their sizes ranging to the 
infinitesimally small. For this special construction, it 
is enough to consider the electric field within a single 
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composite cylinder in a matrix having the effective di- 
electric constant eg = ei[l — 2a?'l3/{a?l3 + 6^)], where 
fi — {€2 ~ ei)/(e2 + ei)- Let the constant applied field 
point in the z-direction, Eo(r) = Eq^- Under this con- 
dition, the presence of the composite cyhnder does not 
change the distribution of the fields in the composite for 
r > b nor the total energy stored in the region occupied 
by the cylinder. 

Within a cylindrical inclusion, the solution of 
with the boundary condition u{r) = —Eorcos{d) reads 



u(r) 



Ar cos(6'), 
(Br+f)cos(0), 



for r < a, 
for a < r < b. 



(4) 



Consequently, the magnitude of the electric field is con- 
stant |E(r)| = \A\ for r < a, and 



|E(r)| 



+ + 2/3(^)2 cos(20) (5) 



for a < r < 6. The coefficients A,B, and C depend on 
the geometry and material properties: 



A = Eoil-p)/if3^2-l), 

B^Eo/i/3^2-l), 

C = Ey 13 1 ((3(1,2 - 1). 



(6) 
(7) 
(8) 



In the study of electric-field fluctuations, it is conve- 
nient to introduce a density of states per unit volume, 
g{E), defined so that g{E)dE is the total number of states 
in the range between E and E + dE, divided by the total 
volume of the inclusion: 



J 5iE~mr)\)dr. 
Substituting into (||) yields 



(9) 



giE) ^ cj)26iE - \A\) 



2(j)2E 



dx- 



1 



where Q{x) is the Heaviside step function, and 

^E{x)^A{(3xr-[l + {(3xf-{E/Br]\ 



(10) 



(11) 



We note that ^e{x) is bounded from above, and that 
the integrand in Eq. ([l0| ) is non-zero only if 7b (x) > 0, 
implying that the density of states g{E) is non-zero only 
for fields in the finite bandwidth: E G [-Ej^^jj^, S'max]- 
The extremal field values for the problem at hand are 
^min = 1^1(1 - and i?max = \B\il + |/3|). 

The highest field determines the electrical breakdown 
properties. Equally important are the fields at which 
g{E) has its maxima, i.e., the fields that occur most fre- 
quently in the composite. To identify them, note that the 
dominant contributions to the integral in (^) come from 
regions close to the points xq = ± |-E/-B|| which 



are solutions of ^e{xo) = 0. Expanding g{E) about a;o, 
and integrating @ yields g{E) ~ (a; - xoY^^H^^^^ or 
g{E) - [(1 - xo)^/^ - {(f)2 - xo)^^^]- Thus, the density of 
states g{E) is not an analytic function since its deriva- 
tives are singular at i? = Evh for which xq ^ 1 or (f>2'- 



E^ 



VH 



^mm' 
|S||1-/302| 

|S||1+/?02| 

^-max- 



(12) 



These singularities are very similar to the van Hove sin- 
gularities found in the density of states for phonons and 
electrons in solids As in the case of solids, the singu- 
larities in g{E) are gernerally associated with saddle and 
inflection points on the surface of generated flelds |E(r)|. 
At these points, g{E) exhibits characteristic sharp local 
maxima or minima, with abrupt changes in the slope. 
This is illustrated in Figs. |]and Because the HS con- 
struction lacks translational symmetry (unlike the subse- 
quent period example), only inflection points are present 
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FIG. 1. The upper half of a composite cylinder of 
outer radius b = 1 with inner core of radius a = 0.5 (i.e., 
4>2 = (ci/b)^ — 0.25), dielectric constant £2 = 10, and outer 
concentric shell of dielectric constant ei = 1: (a) Electric 
field density plot in which lighter shades (gray-scale repre- 
sentation) correspond to increasing field magnitudes, and as- 
sociated contours of equipotential lines u{r) — const, (b) 
Distribution of fields contributing to van Hove singularities 
in the density of field states g{E). Dotted lines show posi- 
tions of Evh = \B\\1 — P(t>2\ ~ Eq while dashed lines identify 
Evh = |-B|(1 -I- Pfpi)- The thick lines denote the boundaries 
of the cylinders. 
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in the density plot of |E(r)| shown in Fig. |l|(a). They 
occur in phase 1 along the lines 9 = and 6 = 7r/2. 

Figure [|(a) shows the density plot of local electric 
fields together with the contours of the equipotential 
lines. Inside r < a region, the magnitude of electric 
field is constant. This generates the i5- function in the 
density of states (|lO|). To locate spatial coordinates 
of the singular fields Evh, we first consider the case 
£2 > ei- From Eq. ^) it follows that Evh ~ -^min 
where {r = a,6 = 7r/2). Since E^^^^^ = \A\, the posi- 
tion of this van Hove singularity overlaps with position 
of the (5-function in g{E). The highest electric fields are 
expected to occur at the two-phase interface. Inserting 
Evh = Euiax into Eq. gives the corresponding coor- 
dinates {r = ,6 = 0) and (r = a+, 6* = tt). Notice that 
at r = a, the surface |E(r)| has finite jump due to the 
discontinuity in normal component of |E| given by (^). 

The locations of the singularities Evh = ± /3(j)2\ 
are not obvious. In this case, the solution of (^) are the 
curves 



6 = ± arccos ( —01- 
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]-(±)^}+fc7r, (13) 



k = 0, ±1,±2,... , and r G [a, 6] which are plotted 
in Fig. |^(b). The dotted central curves correspond to 
Evh = \B\\l ~ l3(j}2\ = E^. The field strength on the 
dashed curves is Evh — 1^1 11 + 134>2\- 

The case in which 62 < ei can be studied using similar 
considerations. For example, the peak of the ^-function 
will be at Evh = -Emax since the highest fields are al- 
ways generated in the less conducting phase. 
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FIG. 2. Density of states g{E) for the HS construction 
with at (1)2 = 0.25 for contrast values a = e2/ei = 2, 10, and 
100. 

In Fig. |[ the density of states g{E) of the composite 
cylinder model is plotted for several values of the contrast 
ratio a = £2/61, at fixed volume fraction (j)2 = (i/^)^ — 
0.25. Increasing the contrast between inclusions leads 
to a broadening of g{E). With an increase of a, the 
maximum field strength (in the less conducting phase) 
increases while in the higher conducting phase, the am- 
plitude of the constant field decreases. The positions of 



van Hove singularities are easily identified as the mini- 
mum or maximum field, or local maxima in g{E). In the 
opposite limit, a ^ 1, the bandwidth of the allowed fields 
collapses to the single peak ^ 5{E — |A|), as expected. 
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FIG. 3. Density of field states g[E) for square arrays of 
cylinders for different contrasts a = ^2/^1, and different vol- 
ume fractions: (j>2 ~ 0.2 (full lines), 0.4 (dotted lines). Poten- 
tial and electric fields are calculated by the finite difference 
method with resolution L/400, L = 1. 

Next we extend our considerations to composite mate- 
rial with periodic microstructure. In particular, we con- 
sider a square array of cylinders of dielectric constant €2 
in a matrix with dielectric constant ei: a model whose 
effective dielectric constant has been well studied [|lO|,0 . 
Let the distance between the centers of neighboring cylin- 
ders be 2L and assume that the external field Eq is ap- 
plied along the x-axis. Although analytical multipole ex- 
pansion techniques (leading to an infinite set of linear 
equations which must be truncated) yield accurate esti- 
mates of the effective dielectric constant, they are not 
adequate to obtain the field distribution. Accordingly, 
we use the finite difference scheme to solve (0)-(||). 

The essence of the numerical method is to map the 
composite to a network of conductors in which each con- 
ducting bond has the dielectric constant of the corre- 
sponding region in the composite. The potential fields 
at each internal node are solutions of the system of 
equation '^jSijiui — Uj) = 0, where index j runs over 
the nodes which are the nearest neighbors of the node 
i, where e^j is the dielectric constant of the bond be- 
tween i^^ and j^^ node. In addition to these equations, 
macroscopic boundary conditions are imposed. Because 
of periodicity, it is enough to consider the solution in a 
square box of size L with a cylinder centered at one of 
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its corners. Then the boundary condition on the edges 
along the direction transverse to the appUed field, say 
the y-direction, is dyu{x,y) = dyu{x,y + L) = 0. In 
the direction of the applied field, there is finite gradient 
u{x + L,y)- u{x, y) = -Eq. 

The solutions for the density of states g{E) are shown 
in Fig. ||. First we notice that in addition to the 6- 
function-like peak associated with the fields inside the 
cylinders, there are three prominent local maxima. These 
maxima are signatures of van Hove singularities. (Fur- 
ther details about the topological properties of the field 



surfaces will be given elsewhere |13 .) The field fluctua- 
tions lie in the bandwidth of the density of states. Fig- 
ure ^ illustrates how by increasing the phase contrast 
a, the fluctuations grow. It seen that the same effect 
is found if the volume fraction 02 increases toward the 
percolation threshold. 

As a flnal example, we consider random arrays of cylin- 
ders whose density of states is distinctly different than 
the first two examples. Previous work on field fluctua- 
tions reveal distributions with two global (as opposed to 
local) peaks 1^,0, although they arose for different rea- 
sons in these two studies. The presence of two global 
peaks is readily understood from our study as is ex- 
plained below. In fact, we conjecture that random mul- 
tiphase composites with widely different dielectric con- 
stants will have many well separated global peaks. 
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FIG. 4. Density of states for a single realization of a ran- 
dom dispersion of nonoverlapping cylinders (see inset) with 
a = 10, 02 ~ 0.4, calculated by a finite difference scheme 
with a resolution L/400. 

Figure ^ shows the histogram of g{E) that we have cal- 
culated for a single realization of a random distribution 
of 20 nonoverlapping disks of dielectric constant €2 = 10 
in a matrix of ci = 1. The system size is L = 1, so 
that the volume fraction (f>2 ~ 0.4. Neglecting the ir- 
regularities due to insufficient statistics in the number 
of calculated fields, the positions of the three local max- 
ima (van Hove singularities) at |E(r)|/|Eo| > 1 are still 
evident, although significantly diminished relative to the 
previous periodic example. Averaging fl^ the density of 
states over many realizations of random samples, as was 



done in Refs. additionally smears the sharpness of 

the band edges and maxima. Therefore, disorder causes 
the local maxima to broaden and eventually merge to- 
gether, i.e., the local features disappear as the disorder 
increases. Similar disorder effects are well known in the 
theory of elecronic density of states in solids. 

To summarize, we find van Hove type singularities in 
the density of states for local electric fields induced by an 
applied field in composites. We show how these singu- 
larities are related to the maxima and the bandwidth of 
g{E). The complex multiple-peak behavior in the density 
of states are not adequately characterized by straight- 
forward calculation of their lower moments. Analysis of 
g{E) at its van Hove singular points represents a powerful 
new approach to quantify field fluctuations in composite 
media. 

It is noteworthy that the density of states analysis of 
field fluctuations laid out in this paper for dielectric com- 
posites can be applied to other field phenomena, includ- 
ing strain fields in elastic media, velocity fields for fiow 
through porous media, and concentration fields for diffu- 
sion and reaction in porous media. 
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